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Abstract 


Algorithms are developed to extract atmospheric boundary layer profiles for turbulence 
kinetic energy (TKE) and energy dissipation rate (EDR), with data from a meteorological 
tower as input. The profiles are based on similarity theory and scalings for the atmospheric 
boundary layer. The calculated profiles of EDR and TKE are required to match the observed 
values at 5 and 40 in. The algorithms are coded for operational use and yield plausible profiles 
over the diurnal variation of the atmospheric boundary layer. 
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Nomenclature 


( = turbulence kinetic energy 

/ = Coriolis parameter 

g = gravitational acceleration 

h = height of atmospheric boundary layer 

L = Obukhov length 

k = von Karman constant, 0.4 

Ri = gradient Richardson number 

7' 0 = reference temperature 

u — mean velocity in the direction of the wind 

v t = friction velocity 

it' , v' = horizontal velocity fluctuations 

w t = conv('ctive velocity scale 

w’ = vertical vtdocity fluctuation 

z = eh'vation above tin' ground 

f = turbulent kinetic energy dissipation rate 

C = Monin-Obukhov stability parameter, zj L 

0 V = mean virtual potential temperature 

0 V ' = virtual potential temperature fluctuation 

4>h - dimensionless potential temperature gradient in the surface lave 

(j) m = dimensionless wind shear in the surface layer 


v 




1 Introduction 


To safely reduce aircraft spacing and increase airport capacity, NASA is developing a 
predictor system, called the Aircraft VOrtex Spacing System (AVOSS: Hinton, 1995; Hinton 
et al.. 1999; Perry et. al., 1997). This system includes prediction algorithms for aircraft 
wake vortex transport and decay. Semi-empirical vortex prediction algorithms have been 
developed and incorporated within the AVOSS (Robins et. al., 1998; Sarpkaya et al.. 2000). 
One of the key input elements for the AVOSS prediction algorithms is the atmospheric 
boundary layer (ABL) turbulence of which the intensity can he represented by turbulent 
kinetic energy (TKE) or eddy energy dissipation rate (EDR). While the prediction algorithms 
require the vertical profiles for tin* TKE and EDR at least up to the vortex generation height, 
the observational data for the TKE and EDR are available only from a meteorological tower 
at the heights of 5 and 40 m above the ground. 

In supporting the NASA AVOSS project, the wake vortex research group at. North Car- 
olina State University (NCSU) has developed algorithms and software which can generate 
the vertical profiles of TKE and EDR. These algorithms are based on the ABL similarity 
relations (Arya, 1988, 1995, 2000; Caughey et ah, 1979; Rao and Nappo. 1998; Sorbjan. 
1989; Stull, 1988) and available experimental data. 

Section 2 describes the ABL similarity relations with respect, to the vertical profiles of 
the TKE and EDR that are dependent, upon the ABL stability. Section 3 contains a detailed 
description of the software. In Section 4, estimates from the similarity relations are compared 
with the observed data. Section 5 provides information on running the software. Finally, 
the summary of this study is given in Section 6. 
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2 Similarity Relations for the TKE and EDR Profiles 


ABL observations frequently show consistent and repeatable characteristics from which 
empirical similarity relationships have been obtained for the variables of interest such as 
TKE and EDR. Similarity theory is based on the organization of variables into dimensionless 
groups that come out of the dimensional analysis. The dimensional analysis is a technique 
used in science and engineering to establish a relationship between different quantities. The 
functional relationships between dimensionless groups are referred to as similarity relations, 
because they express the conditions under which two or more flow regimes would be similar. 
Similarity relationships that have certain universal properties are usually designed to apply 
to equilibrium (steady-state) situations. One of the well-known similarity relations is the 
logarithmic velocity profile law observed in surface or wall layers under neutral stratification. 

The proposed similarity relationships for TKE and EDR are fundamentally based on the 
Monin-Obukhov similarity (Monin and Obukhov, 1954) and mixed-layer similarity (Dear- 
dorff, 1972) theories. The former is applied to a stratified surface layer and is sometimes 
called the surface-layer similarity theory, whereas the latter is applied to mixed layers that 
often develops during daytime convective conditions. These similarity theories have provided 
the most suitable and acceptable framework for organizing and presenting the ABL data, as 
well as for extrapolating and predicting certain ABL information where direct measurements 
of the same are not available. 

Using the framework of these similarity theories, a variety of similarity relationships have 
been suggested to describe the vertical profiles of mean and turbulence fields as functions of 
the dimensionless groups z/L and/or z/h, covering whole ABL including the surface layer 
(Arva, 1988, 1995, 2000; Caughey et ah, 1979; Hogstrom, 1996; Rao and Nappo, 1998; 
Sorbjan, 1989; Stull, 1988). Occasionally, various investigators have suggested different 
values for the empirical coefficients. Based on the similarity scaling in the atmospheric 
surface layer and boundary layer under different stability conditions, the expressions and 
parameterizations for the vertical profiles of TKE and EDR and related characteristic scales 
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are suggested in the following subsections; some of the expressions are adopted directly 
from those references, whereas the others are derived using the similarity relationships of 
turbulence variables other than TKE and EDR. 


2.1 Neutral and Stable Boundary Layers ( zjL > 0) 

The boundary layer may be subdivided into a surface layer (in which stress is nearly 
constant with height) and an outer layer. A separate set of algorithms is assigned to each 
sublayer as follows. 

2.1.1 Surface Layer (z < 0.1//) 


In the surface layer, the TKE (e) and EDR (c) are given bv (Hogstrorn, 1996; Rao and 
Nappo, 1998) 


c = Our. 


f = 


k 


1.24 + 4.3- 


where k — 0.4 is von Karman constant. The friction velocity, //», is defined as 

,21 1/2 


lit = [(«'«•')? + («'«■’')* 


( 1 ) 

(2) 

(3) 


where the right hand side of Eq. (3) represents the total vertical momentum flux near the 
surface (the subscript. ,s denotes the ground surface). The Obukhov length L depends on 
both the momentum and heat fluxes near the surface and is defined later: the ratio zjL is 
the fundamental similarity parameter of the Monin-Obukhov similarity theory. 


2.1.2 Outer Layer (c > 0.1/?) 


Expressions for the outer layer can be assigned according to the level of stratification. 
(1) Neutral and Stable Boundary Layer 

In the neutral and moderately stable boundary layer, the TKE and EDR are given by 
(Hogstrorn, 1996; Rao and Nappo, 1998) 

—•Hr 
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( 5 ) 


E “l( 1 ' 2J+4 ' 3 i)( 1_a85 i) - 

Alternatively, Eqs.( 4) and (5) may also be used for the entire boundary layer, including the 
surface 1 layer. 

(2) Very Stable and Decoupled Layers 

In the very stable boundary layer and decoupled layers, the TKE and EDR can be 
expressed by extension of Eqs.(l) and (2) as 

e = 6 u\, ( 6 ) 


where u L is the local (friction) velocity scale and L l is the local buoyancy length scale. Under 
very stable conditions, the elevated layers of turbulence are decoupled from the surface and 
the local fluxes and scales cannot be reliably estimated. Perhaps, an empirical relationship 
between the overall turbulence intensity (e l /‘ 2 /u) and Richardson number should be explored. 
It is worthwhile to note that some experimental results show that Eq.{ 5) can be still used 
to estimate f even in a very stable boundary layer. 

2.2 Unstable Boundary Layer (z/L < 0) 

The unstable ABL such as during daytime surface heating can be divided into three 
regimes depending upon the stability parameter, z/L or h/L. 

2.2.1 Strongly Unstable (Convective) Regime (\z/L\ > 0.5) 

The structure of the convective regime is dominated by buoyancy. The mean wind 
velocity and potential temperature profiles are nearly uniform with height. For this reason, 
the convective outer layer is called the “mixed layer.” The mixed layer is topped by an 
inversion layer in which temperature increases with height. A broad maximum of TKE is 
usually found in the middle of the mixed layer, while EDR decreases slightly with height. 
(1) Surface Layer (2 < 0.1/*.) 
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( 8 ) 


In the surface layer, the TKE and EDR are given hv (Area. 2000) 

e = 0.36<c* + 0.85«* fl — 3— 
u;V 2 2/ 3 Y V2 

, = *i( 1+a5 z ) 

(2) Mixed Layer 

In the mixed layer, the TKE is given by (Arya, 2000) 


r = [ 0.36 + 0.91 f] ( 1 - 0.8 
Jr 


)d' 


or, for most practical purposes, 

r = 0.54 uJ. 

The EDR decreases slowly with height at a linear rate (Sorbjan, 1989), i.e.. 


r =f r-° 

where 4 the convective velocitv scale w * is defined as 




a \ ^ 

J() / 


(9) 


Here g is the gravitational acceleration, T 0 is the reference temperature, and (w'0 v ') s 
mean surface heat flux. 


( 10 ) 

(11) 

( 12 ) 

(13) 

is t.h(* 


2.2.2 Moderately Unstable Regime (0.02 < \z/L\ < 0.5) 

In this regime, the mechanical production of TKE is comparable with buoyancy pro- 
duction of TKE, i.e., turbulence generation from vertical wind shear is comparable to that 
generated from surface heating. The TKE in this regime is more or less uniform over the 
boundary layer or may decrease slightly with height, and the boundary layer structure may 
be more uncertain. 

2.2.3 Weakly Unstable (Near-Neutral) Regime (\~-/L\ < 0.02 or |/t/L| < 1.5) 

This regime often occurs during the transition period of early morning and late afternoon 
or during overcast days with strong winds. The lapse rate for temperature tends to be near- 
neutral. In this regime, mechanical (shear) production dominates the TKE budget. 



3 Software Description 

In order to generate vertical profiles of TKE and EDR for operational applications, 
software is written which utilizes the algorithms in the previous section with observations 
measured at 5 and 40m above the ground. With the determination of the characteristic 
similarity scales (such as L, «*, w„ and h), the computation of the TKE and EDR profiles 
is straightforward. The characteristic similarity scales can be estimated from the winds and 
virtual potential temperatures measured at two levels near the ground. 

Since the atmospheric similarity relationships are based on the mean quantities of the 
observed wind, temperature, and turbulence, the required time averaging interval for EDR 
and TKE profiles should be at least 30 minutes (Stull, 1988). Hence, with the exception of 
TKE, all measured variables are 30-minute averaged. For the measured TKE, a 30-minute 
median value is used, since sporadic measurements of exceptionally large TKE can cause an 
unrealistic TKE average. 

The vertical profiles generated from similarity theory are additionally required to match 
the measured values at z = 5 and 40 m as closely as possible. Often the profiles generated 
from similarity expressions do not exactly match the observed values at z = 5 and 40 m 
simultaneously. Therefore, the TKE and EDR profiles between 2 = 5 and 40 m are assumed 
to be linear and independent, of the similarity relationships. In addition, upper and lower 
bounds for the values of TKE and EDR at , z = 5 m are specified to prevent unrealistically 
large or small values compared with those at z = 40 in. Above 40 m, the similarity profiles 
are used, but are adjusted to match the observed TKE and EDR at 2 = 40 m. 

3.1 Determination of Surface Layer Similarity Scales 

According to algorithms in the previous section, software requires to first determine 
surface layer similarity scales L and «». 

The Obukhov length is defined as 

L = - W- , (14) 

k(g/To)(w%X 

I 
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in which friction velocity and the surface heat flux can be estimated from measurements of 
tin' mean differences or gradients of velocity and temperature between any two heights cj 
and z -2 within the surface layer, but well above the tops of roughness elements. 

Letting Au = « 2 — «i and Al9„ = 0 V 2 — 0 V 1 be the difference in mean velocities and 
virtual potential temperatures across the height interval A; = z 2 — Cj , one can determine 
the gradient Richardson number (Ri) at the geometric height z m = (~i~ 2 ) I/2 by 


RH 


y I , -2 

= —Zr„ In — 


A 0 V 


\2 ' 


To V zj (Any 
The corresponding value of the Monin-Obukhov stability parameter Q, 
determined from the relations given by (Arya, 1988) 


(15) 


,[L can be 


Cm = Ri{z m )i for Ri < 0 , 


(1G) 


Cm 


for 0 < Ri < 0.2. 


1 - 5 Ri(z m ) 

Then, the basic universal similarity functions (j) ni and cl>i, are directly related to Cm-, i e., 


(17 


<t>h = <P,n = (1 - 15Cm) 1/2 • for Cm < 0 O 8 ) 

(f>h = 4>m = (1 + 5Cm). for Cm >0 (19) 


The similarity relations of Eqs. (17) and (19) are not valid for Ri > 0.2 and Cm > 1- 

Finally, the friction velocity and the surface heat flux can be obtained from the following 
relations 


= 


k Au 

®m(Cm)ln(z 2 / Zi)' 

k 2 AuAffr 

Vm(Cm)$h (Cin) (jn-^^ | 


( 20 ) 

( 21 ) 


The parameters Ri(z m ) and z m /L are used to determine the stability criteria discussed 
in the previous section, while the friction velocity u» is used to estimate the ABL height h 
in the neutral and stable boundary layer. 



3.2 TKE and EDR Profiles in Neutral and Stable Boundary Layer 

(z m /L > 0) 


The boundary layer height can be estimated as the minimum of those given by the 
following diagnostic relations (Arya, 1995): 


h = 0.3 


u * 


r 


h = 0.4 1 


u„L 

~T. 


1/2 


( 22 ) 


(23) 


where / is the Coriolis parameter. Note that Eq.( 22) is valid only for a stationary neutral 
boundary layer, but it gives an upper bound for h in slightly stable or near-neutral conditions 
when L becomes too large and Eq.( 23) overestimates h. 

Then, Eqs.( 4) and (5) are applied to obtain TKE and EDR profiles above 2 = 40 m 
height. The friction velocity in Eq. s.(4) and (5) is adjusted so that the profiles are continuous 
at z = 40 m, and might be slightly different from that given by Eq. (20). 


3.3 TKE and EDR Profiles in Unstable Boundary Layer (z m /L < 0) 

3.3.1 Strongly Unstable (Convective) Regime (\z m /L\ > 0.5) 

In this regime, the ABL height is usually estimated from the height of the inversion base, 
which can be determined from the vertical temperature sounding. However, it is difficult to 
obtain the temperature inversion base from the existing sounding data from the Dallas/Ft. 
Worth (DFW) and Memphis field experiments. This is because the soundings often do not 
extend above a height of 1 Arm, whereas the ABL height in this regime can easily reach a 
height of 2 — 3 km. Alternatively, we estimate the boundary layer height from £>/.*>■. (1 1 ) and 
(12), using the measured TKE and EDR at z = 40 m, i.e., 

0.4 -I- ^0.16 + 0.3 240 f4o/( e 4o/0.54) 10 

h = 77 7777770 (24) 

Uo/ (£40/0.04) 

where subscript 40 represents 2 = 40 m. Although the height of z — 40 m is within the 
surface layer (z < 0.1 h) in most of cases of this regime, the observations indicate that the 
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mixed layer similarity may also he used for the entire boundary layer except for the layer 
very close to the surface. 

Then. £</.s.(10) and (12) are applied to obtain TKE and EDR profiles above c = 40 m 
height., respectively. The convective velocity scale in £7/. (10) is adjusted for the TKE profile 
to be continuous at z = 40 m, and will be slightly different from that estimated from £7/. (11). 

3.3.2 Moderately Unstable Regime (0.02 < \z m /L\ < 0.5) 

This regime also includes the cases of c r) > <^40 but \z rn /L\ > 0.5, since TKE (i.e.. c) 
usually increases with height near the surface in the convective regime. The profiles are 
obtained in the same wav as in the convective regime except that £r/.(ll) is used for TKE 
profile rather than £</.(10). 

3.3.3 Weakly Unstable (Near-Neutral) Regime (\z m /L\ < 0.02 or \h/L\ < 1.5) 

I11 this regime, the formulations for the neutral boundary layer are applied for the profiles. 
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4 Comparisons with Observed Data 


The friction velocity and TKE and EDR at 2 — 5 m and 40 m estimated from the simi- 
larity relations in the previous sections have been compared to measurement data from the 
Dallas/Ft. Worth (DFW) Airport Field Experiment (Dasey et al., 1998) during the period 
of September 15 - October 3, 1997. The height of 2 = 5m is essentially assumed to be within 
the surface layer. The measured u* and EDR are 30-minute averaged, while the measured 
TKE is a median value over 30-minutes. The frequency for each stability regime is shown in 
Table 1, indicating that most of data belong to stable and moderately unstable regimes. In 
the following, the comparisons are performed in two stability groups, that is, neutral and sta- 
ble regime including the weakly unstable; regime and unstable regime including moderately 
unstable and convective regimes. 

4.1 Friction Velocity (u*) 

Figure 1 shows that the estimated w* from mean velocity measurements with Eq. (20) 
agrees well with the measured u* from the momentum flux at 2 = 5 m, but the data scatter 
becomes larger for weaker it*. In particular, for very small values of it*, which are often 
obtained during very stable conditions, the estimated u t tends to be significantly smaller 
than the measured u». This is because the Monin-Obukhov similarity theory on which it, 
estimates are based, is not applicable under very stable conditions (Nappo and Johansson, 
1999). 

4.2 TKE at z = bm and 40 m 

As formulated in sections 2 and 3, the estimated TKE at 2 = 5 rn in neutral and stable 
conditions is given as a function of only it,, while in unstable conditions, it is given as a 
function of it*, w t , L and 2 . As shown in Fig. 2, the estimated TKE appears to agree with 
the measured, although data scatter becomes larger for smaller values. This is especially 
true in neutral and stable conditions. For very small values of TKE, most of the estimated 
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TKE values are significantly lower than the measured TKE. The agreement between the 
estimated and measured TKE is much better for the unstable conditions. 

At 2 = 40 m, only a comparison for the neutral and stable conditions is conducted. 
In unstable conditions, the ?/;, and h are calculated from the measured EDR and TKE at 
2 = 40 m and thus, the estimated TKE and EDR at 2 = 40 m will match their measured 
values at 2 = 40 in. exactly. Figure 3 shows that similar to TKE at 2 = 5 m, the estimated 
TKE in neutral and stable conditions agrees reasonably well with the measured TKE, but 
data scatter at 2 = 40 m is larger than that at 2 = 5 m, especially for smaller TKE values. 

4.3 EDR at z = om and 40m 

Using the estimated values of u* and L. EDR at 2 = 5 m is calculated from Eqs. (2) and 
(9). Its ‘measured’ values were obtained from the power spectrum of the observed wind data 
using the theoretical Kolmogorov relations of the spectra in the inertial subrange. 

Figure 4 shows that at 2 = 5 m, the estimated EDR appears t.o increase linearly with 
increasing measured EDR. but is considerably overestimated in both stability regimes. Data 
scatter is somewhat larger compared to that of TKE. Similar to that of TKE. however, the 
data scatter increase's with decreasing values of EDR and is larger in neutral and stable 
conditions than in unstable conditions. As explained in section 4.2, the comparison for EDR 
at 2 = 40m is conducted only for the neutral and stable conditions. Figure 5 shows that 
agreement and data scatter between estimated and measured EDR at 2 = 40 m are similar 
to those at 2 = 5 m, but with smaller overestimation. 

4.4 Adjustments in TKE and EDR Profiles 

As shown in the above comparison plots, the ABL similarity theory appears to represent 
the measured 1 /* and TKE at 2 - 5 m and 2 - 40 m reasonably well, but with increasing data 
scatter for decreasing values of and TKE. On the other hand, the estimated EDR consid- 
erably overestimates the measured EDR, although estimated and measured values appear to 
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he strongly correlated with a linear relationship between the two. As described in section 3, 
therefore, we require that the generated profiles for TKE and EDR match the measured val- 
ues at 2 = 5 and 40 m while maintaining the profile shape from the similarity theory above 
z = 40 m. Otherwise, the estimated profiles for TKE and EDR may significantly deviate 
from those ‘measured’, especially for small values of TKE and EDR. 
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5 Procedure of Running Software 


The software is coded in FORTR AN and is designed for operational use. Tlx' code is 
composed of one main program and two subroutines. As input the main program reads data 
files containing TKE and EDR at 2 = 5 m and 2 = 40 m, as well as mean winds and virtual 
potential temperatures at two levels near the surface. The subroutine “EDRTKEJSTABLE” 
computes TKE and EDR profiles under neutral and stable conditions. The other subroutine 
“EDRTKE_UNSTABLE” computes TKE and EDR profiles under unstable conditions. O 11 
output, vertical profiles of TKE and EDR are generated from the main program. 

In executing the code, one needs to enter the following: 

• latitude (degree) of the airport , which is used to calculate the Coriolis parameter /. 

• one-minute averaged tower data file for wind and virtual potential temperature at two 
levels near the surface. 

• measured 30-minute averaged EDR data file at 2 = 5 rn. 

• measured 30-minute averaged EDR data file at 2 = 40 in. 

• measured 5-minute averaged TKE data file at 2 = 0 in. 

• measured 5-minute averaged TKE data file at 2 = 40 in. 

Once the input data files have been successfully entered, the code then generates an 
output file, which contains not only the vertical profiles of TKE and EDR . but also the 
output of similarity scales such as h. L. and ?/?*. Figures G and 7 show a diurnal variation 
of TKE and EDR profiles generated by the software for a typical sunny day. 
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6 Summary and Discussion 


Based on the existing similarity theories in the atmospheric surface layer and bound- 
ary layer under different stability conditions, the expressions and parameterizations for the 
vertical profiles of TKE and EDR have been suggested. Compared with observation, theoret- 
ically estimated TKE at 5 and 40 m heights above ground appears to agree reasonably well 
with those observed at the same heights except for the very small values of TKE. However, 
theoretically estimated EDR at 5 and 40 m heights considerably overestimates the observed 
EDR with large data scatter, although estimated and measured values appear to be strongly 
correlated with a linear relationship between the two. 

From the ABL similarity relationships and their comparisons with the observations, soft- 
ware has been developed to generate realistic vertical profiles of TKE and EDR. The input 
parameters for the software are the measured TKE and EDR at the heights of 2 = 5 and 
40 rn above the ground, and the measured winds and virtual potential temperatures at two 
levels near the ground from which characteristic similarity scales, such as L and u*, can be 
estimated. In the software, to minimize the difference between the similarity relations and 
observations it has been required that calculated values match the observed values at the 
heights of 5 and 40 m while maintaining the profile shape from the similarity theory above 
40 m. Although the software yields very plausible vertical profiles and their diurnal varia- 
tions, statistics for the difference between the values estimated from the software and the 
measured values of TKE and EDR at the heights other than 5 and 40 m has to be obtained 
for further improvement of the software. 
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Stability 

Frequency (%) 

Stable 

54 

Weakly Unstable 

5 

Moderately Unstable 

37 

Strongly Unstable 

4 


Table 1: Stability frequency (%) in the Dallas/Ft. Worth (DFW) Airport Field Experiment 
during the period of September 15 - October 3, 1997, where data from two rainy days and 
one day for which data are partly missing have been omitted from the total data set. 
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